Effect of Time Reversal Symmetry Breaking on the Density of States in Small 

Superconducting Grains 
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We show that in ultra-small superconducting grains any concentration of magnetic impurities or 
infinitely small orbital effect of magnetic field leads to destruction of the hard gap in the tunneling 
density of states. Instead, though exponentially suppressed at low energies, the tunneling density 
of states exhibits the "soft gap" behavior, vanishing linearly with excitation energy, as the energy 
approaches zero. 

PACS numbers: 73.23.b, 74.80. Bj 



The gap in the tunneling density of states (DoS) is one 
of the most fundamental manifestations of microscopic 
mechanism Q| of superconductivity ^]. The existence 
of this gap is closely related to the time-reversal sym- 
metry in the superconductor. As a result, non-magnetic 
impurities (at low densities) do not affect the gap, A, 
in conventional superconductors (the Anderson theorem 
||). It is magnetic impurities that violate the symmetry 
of the superconducting state and act as pair breakers. 
Therefore bound states below the gap are created and 
the gap is suppressed Q. This effect of the suppression 
of the gap has been observed experimentally and it is 
now a textbook example || . 

If the density of the magnetic impurities n m exceeds 
certain critical value, n c the mean field transition be- 
comes forbidden at all. Even at much smaller density, 
the magnetic impurities suppress the DoS. This suppres- 
sion of the hard gap manifests itself in the appearance of 
the finite DoS at energies (measured from the Fermi level) 
e, smaller than the mean field gap |e| < A, [Q. However, 
at n m < 0.9n c , the self-consistent Born approximation 
(SCBA) of Ref. Q predicts the DoS still vanishing at en- 
ergies smaller than a certain value, which is referred to 
as the renormalizcd gap e*, (see dashed line in Fig. 1.) 
In the region 0.9n c < n m < n c , the so-called regime of 
gapless superconductivity occurs, where the gap in DoS 
vanishes even though the superfluid density still remains 
finite. 

The prediction of SCBA about the existence of the 
finite gap in the density of states even at the already 
broken time reversal symmetry is intriguing since there 
is no longer a symmetry reason for vanishing of the DoS. 
The situation in some sense reminds the problem of ex- 
ponentially small tails in the DoS at small energies in 
normal materials The appearance of such tails is 

non-perturbative. It means, that any approach which 
consists of selecting some class of diagrams based on ex- 
pansion in disorder strength will result in zero DoS under 
the gap. Rather, in order to find the shape of such tails 



one has to perform the instanton analysis (also known as 
the "optimal fluctuation method" [Q-[|). In this Letter, 
we perform such analysis for the effect of magnetic im- 
purities or small magnetic field on the DoS of supercon- 
ducting grains |5|. We will demonstrate that the DoS at 
< |e| < A becomes finite, though exponentially small, 
no matter how small the magnetic field or the concen- 
tration of magnetic impurities is. At |e| — > 0, the DoS is 
shown to exhibit the "soft gap" behavior vanishing lin- 
early with excitation energy. 

The Bogolyubov equations |To| for quasiparticle spec- 
trum in s-wave superconductors are 



eip = Hip, 

where tp is the Gor'kov-Nambu spinor jul 



u a (n) 



(1) 



(2) 



[Unless stated otherwise, Latin (Greek) indices label the 
orbital (spin) states] . The mean field Hamiltonian is 
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where the one-particle Hamiltonian H may act both on 
the orbital and spin coordinates of the electron, and T is 
the time inversion operator, Tu a {n) — c^u^n), where 

a v is the Pauli matrix, T 2 = 1. 

The crucial characteristic of the system is the symme- 
try of H with respect to the time inversion, 

H = H S + H a ; TH S T = H s : TH a T = -H a . (4) 

If the time inversion symmetry is preserved, H a = 0, 
both diagonal entrees of the Hamiltonian (||) can be di- 
agonalized simultaneously, 

iT =e = diag(a,6,---), (5) 
and one obtains the eigenvalues of Hamiltonian H s 



I 



(6) 



so that the hard gap in the one-particle excitation spec- 
trum exists independently of further model assumptions 
on H s J3|. However, if the time inversion symmetry is 
broken, the answer is not universal and we need to fur- 
ther specify the model. We adopt H s and H a to be 
independent M x M, (M — > oo) Random Matrices |12| 
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satisfying the constraint (Q). (We will omit the spin 
indices where it does not cause any confusion.) Here 
t5i = — £j) <C A is the mean level spacing, the 

parameter r^ 1 < A characterizes the strength of pair- 
breaking potential (see below), and (. . .) stands for the 
ensemble averaging. Therefore, H s belongs to either or- 
thogonal or symplectic ensembles and H a describes the 
crossover to the unitary ensemble. In what follows, we 
choose the basis of the eigenstates of the Hamiltonian 
H s , so that H s has the form (||) while H a in this basis 
is a random matrix with the correlation function ((?]) . 

The DoS in the system is expressed in terms of the 
disorder averaged Green function 

(v(e)) = -±-ImTv{G R (e)), G R (e) = } (8) 
2vr e - U + lO 

If the time reversal symmetry is preserved, th -> oo, the 
DoS is given by the usual BCS expression 

Before we proceed, let us discuss the physical situa- 
tions for which the random matrix description (RMT) (Q) 
is applicable. It is well known that RMT description 
of the spectrum requires all the relevant energy scales to 
be much smaller than the Thouless energy Et- It trans- 
lates to the condition that the size of the grains R to be 
much smaller than the coherence length £, see Ref. Q 
for recent experiments on such grains. The parameter 
t^ 1 , can be related to the physical characteristics of such 
grains. If magnetic field is applied it penetrates through 
the small grain without screening and 

t h 1 = E T ($/$ ) 2 , 

where <J> is the magnetic flux penetrating through the 
grain and <&o = hc/2e is the flux quantum. In the case 
of the doping by magnetic impurities th is given by the 
spin-spin scattering time 

= t- 1 = ^ Q N S \V S \ 2 S(S + 1), 

where n s is the concentration of magnetic impurities, S 
is the impurity spin, V s is the scattering matrix element 



and Vq is the thermodynamic density of states. Finally, 
the possibility to neglect the non-Gaussian correlations 
of the Hamiltonian H a is guarded by the requirement 
that no matrix element of H a exceeds the Thouless en- 
ergy. As we shall see below, the characteristic value of 
the matrix elements contributing to the low energy tail 
of the DoS is of the order of A <C Et, so that the Gaus- 
sian approximation is justified. (All of these assumptions 
definitely break down in bulk systems R ^S> £.) 

The exponentially small tails in the DoS were consid- 
ered in Refs. The idea is to look for the fluctua- 
tions of the random potential H a which form a low en- 
ergy bound state, thus leading to non-zero DoS at such 
energy. The probability to form such a bound state is de- 
termined by the distribution of the matrix elements HQ 
and, while exponentially small, should be maximized by 
choosing the "optimal" fluctuation of the potential. The 
resulting DoS is then proportional to this probability. 

To gain some intuition about the form of the optimal 
fluctuations, let us consider the simplest possible realiza- 
tion of the random potential H a where it couples only 
two eigenstates iq, jo of the Hamiltonian H s : 



iV; i = i , j =jo; 
HQ = { -iV; i = ja, j = «o; 

0: otherwise. 



(9) 



All states with i ^= io,jo decouple, and the relevant 
Hamiltonian reduces to 4 x 4 matrix: 

^(0) _ ( 1~Lioia H-ioja 

Ui j =iVlN, Hj i =H* ojo (10a) 
n i0H> = & af + Aaf ; H jojo = + (10b) 

where In is the 2x2 unit matrix in Gor'kov-Nambu 
space, and a^ z are the Pauli matrices in this space. The 
eigenvalues of the matrix ([h]) are 



£1,2 = ± 



A 



(ii) 



e 3A = ±[A+\V\+0{Zl)], £, 



Two of these eigenvalues, £3,4, lie above the gap, and are 
not interesting for our purposes. The other two, £1,2, 
correspond to the bound states under the gap. Clearly, 
by suitable choice of V (and £'s), we can tune ei.2 to 
the desired energy e. The averaged DoS is proportional 
to the probability P(V) to find such value of the matrix 
element. According to Eqs. (Q) and (^|), we have 

P(V) oc exp (-t^ 2 

The minimal value of V, providing the level ( |Tl| ) to have 
energy e, is |V| = A — |e|, and we obtain 



2 



(Ke)>«exp -?(|e|-A) 



(12) 



where we omit all pre-exponential factors (to be calcu- 
lated below). 

Equation ( |l2| ) is the main physical result of this Letter. 
We have shown, that the DoS in the ultra small super- 
conducting grain possesses the exponentially small tail at 
low energies even below the renormalized value of the gap 
obtained using the self-consistent Born approximation. 

To make our derivation rigorous we have (i) to prove 
that the ansatz (^) is indeed a saddle point in the ensem- 
ble averaging; (ii) to calculate the pre-exponential factor 
by summing over all saddle points (io, jo) and integrating 
over the fluctuations around the saddle point. 

Saddle point — To find the saddle point one has to 
minimize the exponent of the Gaussian probability (Q) 



ly 

2 ^ 



(13) 



with respect_to all the matrix elements H^, subjected to 



constraint 



and the condition 



(14) 



where eo{H a } is the smallest eigenvalue of the Hamilto- 
nian (111) . This involves finding a solution of the equations 







i'<j> 



= 



(15) 



for all 1 < i < j < M. Here hyy and A are the Lagrange 
multipliers to be found from the conditions (0) and (14). 
Excluding A^, we find from Eq. (|l|) 



zja 



de <9e 



(16) 



where 



constitute the matrix analogue of the non-linear 
Schrodinger equation of Ref. || . 

The essential simplicity of the random matrix model 
( |l7j ) stems from the independence of the non-linear terms 
A, B of the state index j. Namely, that Eq. ( |l6| ) can be 
considered as a linear equation for a state j, while the 
coefficients A, B have to be found self-consistently. With 
non-linear terms A, B fixed, the non-trivial solution to 
Eq. (H|) for a given eigenvalue e exists only for two val- 
ues of £j [similar to Eq. @]. It means, that at most 
only two states can be mixed. According to Eq. Jl6|), 
it indicates that ansatz ([io]) we adopted from the very 
beginning is the only possible form of the saddle point. 

Pre-exponential factor — Having convinced ourselves, 
that we have found the optimal fluctuation and, thus, the 
exponent (jl^) correctly, we turn to the calculation of the 
pre-exponential factor in this expression. Our starting 
point is once again the optimal fluctuation ([l0|), how- 
ever, we wish to take into account all the other matrix 
elements which couple the states m = io, jo with all 
other states (j ^ io, jo)- Since the exponential tail comes 
from the two levels io, jo, coupled by the large matrix 
element V » 1/^/7, the elements may take only its 
typical value ~ l/\/7 ^ A. Therefore, all of them can 
be treated in perturbation theory. To the second order 
in H^j the effective Hamiltonian acting in the reduced 
Hilbert space of the two states io, jo acquires the form 
dio| ) with the entrees changed due to the mixture of all 
the other levels (indices m,n = io, jo) '■ 
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(18) 



Substituting Eq. (|18|) into Eq. (|3|) we obtain the contri- 
bution of two lowest levels into the (non- averaged) DoS: 



v(e) = \e\(l-Z)6 \e\\-Zf 



W 



(19) 



W 



(1 + Z)A+I±-\H9. 



e+Y 2 



de 



u*(i)u(j) + v{i)*v(j). 



Here u(n), v(ri) are the components of the Nambu spinor 
ipj [see Eq. (||)], corresponding to the eigenstate eo- Sub- 
stituting Eqs. © and (|) into Eqs. (§]), @, we obtain 



^a z +Aa x +A ) tpj - Bi(j 



n n 

which, together with the normalization condition, 

E^ = 1 



(17) 



where £± are given by Eq. (|11[), and we introduced 

M "" j "j'-j 
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For the optimal fluctuation (|J), the argument of the S 
-function reproduces the spectrum (J5~T|) . 

The DoS (|l9|) should be averaged over the fluctuations 
of the matrix elements and summed over states io, jo 



<^)> = e n 



»o Jo i<3 L' 



exp 



1/(6). 

(20) 
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At energies e > ii, we can neglect the level repulsion 
and replace the sum over io,Jo by the integral 



E 

»o Jo 



2^2 / ^£»o^£jo> 



where factor of 1/2 excludes double counting the same 
configurations (io,jo) and (Jo,io). Then straightforward 
integration in Eq. ( p0| ) utilizing the condition |e| 
5i/(Ath) yields the averaged DoS: 



(Ke)) = -^i(e)exp 



(21) 



Fi(e) 



2|e| 



*2(e) = rjf 



ft (A + 4|e|r H VA 2 -e 2 ) 
(A-|£|) 2 
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A - 



Si 




Equation ( pl| ) is the main quantitative result of this 
Letter. It gives the paramctrically exact description of 
the exponential tail in the DoS. It is valid as long as the 
exponent Fi is larger than unity. At low energies e <C A 
we can neglect the second term in Fi(z), reproducing the 
qualitative result Eq. Jl2]). 



v(e)8! 




FIG. 1. DoS for a superconducting grain with broken time 
reversal symmetry. The thick solid line is the tail of the 
DoS (J2 1|) . The dashed line represents the result of the 
self-consistent Born approximation Ref. The inset shows 
the behavior of the tail at e — > 0. The curves are plotted for 
A = 3.55i and t h = 1.557/ 1 . 

At larger energies, when e approaches A, the valid- 
ity of our considerations and of Eq. (|2l]) breaks down at 
the point e* where the two terms in the exponent ^(e) 
become of the same order. Remarkably, that the point 
A— |e*| ~ A/(Arir) 2 / 3 ] is parametrically the same as the 
"renormalized gap" , e* , predicted by SCB A H . Combin- 
ing these results, we can describe the DoS in the grain 
by the continuous function depicted on Fig. 1. 



At small energies e — > result (|21|) vanishes linearly 
due to the pre-exponential factor Fi(e). The appearance 
of this linear suppression is easy to reveal already on the 
level of qualitative analysis. Indeed, one can see from 
Eq. (|l^) that the levels near the e = repel each other 
due to the difference in energies of the one electron state 
£. It means that the contributions of the levels is always 
limited by |£, — £j| < |e| which gives the corresponding 
smallness in the integration domain. Finally, |e| < Si, the 
level repulsion |l2j between those orbits should be taken 
into account which results in the additional suppression 
by a factor (e/Si)" , where (3 = 1(4) in the absence (pres- 
ence) of the spin-orbit coupling. 

To conclude, we have shown that in a small supercon- 
ducting grain breaking the time reversal symmetry leads 
to the appearance of the exponentially small tail in the 
DoS at energies smaller than the BCS gap. The DoS 
is non zero for all energies |e| < A except for the point 
e = 0, where the DoS vanishes linearly in energy. Thus 
the grain no longer exhibits the hard gap in the excita- 
tion spectrum as predicted by Ref. Q]. Our results are 
non perturbative as they were obtained by considering 
the optimal fluctuation of the random potential. 
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